Nodal-antinodal dichotomy from pairing disorder in d-wave superconductors 
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We study the basic features of the local density of states (LDOS) observed in STM experiments 
on high-T c d-wave superconductors in the context of a minimal model of a d-wave superconductor 
which has weakly modulated off-diagonal disorder. We show that the low and high energy features of 
the LDOS are consistent with the observed experimental patterns and in particular, the anisotropic 
local domain features at high energies. At low energies, we obtain not only the scattering peaks 
predicted by the octet model pj, but also weak features that should be experimentally accessible. 
Finally, we show that the emerging features of the LDOS lose their correspondence with the features 
of the imposed disorder, as its complexity increases spatially. 



In recent years, the effect of underlying inhomo- 
geneities in superconductors has been studied inten- 
sively. In the context of high temperature supercon- 
ductors (HTSC), checkerboard local density of states 
(LDOS) oscillations and strong nanoscale gap inhomo- 
geneity have been observed in scanning tunneling spec- 
troscopy (STS) experiments [2], and signals in dynamical 
susceptibility measured by neutron scattering have been 
interpreted as stripe-like nanoscale modulations of charge 
and spin degrees of freedom [3[ . Even though it is still not 
clear whether these modulations are intrinsic or driven 
entirely by disorder, it seems plausible that the magni- 
tude of T c might be strongly related to the very existence 
of inhomogeneity on the coherence length scale 0, [Hj], 
and therefore a deeper understanding of the source of 
inhomogeneities is important. In the cases of the com- 
pounds La2- a 3a ;E Cu04and La2- cc Sr ;E Cu04(with Nd or 
Eu co-doping), inhomogeneities take the form of static 
long-range stripe-like spin and charge modulations • 

Tunneling spectroscopy has been used to probe states 
in different regions of momentum and energy, by the 
help of the Fourier transform scanning tunneling spec- 
troscopy (FT- STS) and the high energy features of the 
LDOS. At low energies, in the d-wave superconduct- 
ing state near optimal doping, quasiparticle interference 
(QPI) observed by FT-STS is dominated by peaks at 
well-defined wavevectors which are consistent with 
a simple model of Bogoliubov quasiparticles, called the 
octet model ToL Hl|. Even though the qualitative fea- 
tures of the octet model are experimentally robust, a 
quantitative understanding of the amplitude, location 
and width of the peaks is not straightforward to obtain 
and depends rather sensitively on the nature of the scat- 
tering medium [l2j]. A point-like scatterer, for example, 
in an otherwise homogeneous dSC leads to a landscape of 
some spot-like and some arc-likedispersive features close 
to in the FT-STS images [13], whereas experimen- 
tal data appear mostly spot-like. At high energies, local 
unidirectionality has been observed in domains with size 
approximately 5ao, close to the superconductor's coher- 
ence length and the domains typically alternate in ori- 



entation. This behavior is usually attributed to a disor- 
dered charge-density wave (CDW), with some success on 
demonstrating the emergence of the LDOS patterns from 
the underl ying CDW order. [141 ]. 

In Refs. |l5l flf^ it was shown that generic off-diagonal 
disorder (coupling of the nodal d-wave quasiparticles to 
an s-wave order parameter A s ) is a relevant operator 
that leads to a first-order transition in a system with 
tetragonal-orthorhombic symmetry breaking. Therefore, 
off-diagonal disorder presents the dominating effect for 
superconductors with tendency towards orthorhombic 
symmetry breaking. However, when the orthorhombic 
symmetry is not broken (due to frustration effects), it 
is expected that modulated forms of the expected dis- 
order should be present, with no effect at large length 
scales. In this paper we present a an explicit model with 
off-diagonal, modulated disorder, which unifies the ba- 
sic features of the STM observations into a consistent 
framework. 

In this model, itinerant fermions are coupled to weak, 
off-diagonal disorder, modulated in well-defined domains. 
Such a state is characterized by a density wave or mod- 
ulated nematic order. We motivate the existence of such 
disorder on the basis of: i) a phenomenological theory of 
competing s and d-order parameters ii) recent mean-field 
studies on more detailed states, supporting a modulation 
of the off-diagonal disorder [14| . First, we show by us- 
ing first order perturbation theory, that the spectrum at 
low energies and the position of the nodes remains un- 
changed, an expected result since at long length scales 
the disorder averages to zero. Secondly, we find that at 
low energies the Fourier-maps of the LDOS have peaks 
that are dispersing according to the octet-model scat- 
tering [IH, given that A^/A d <C 1 with the addition 
of weak static peaks, bearing strong similarities to the 
STM experiments' observations In addition there are 
static peaks close to (±37r/5, ±7r/5), but these peaks are 
very weak (since A^ s ^ is very small) and if spatial disor- 
der exists, span a large region around these wavevectors. 
Thirdly, we show that at high energies the LDOS ac- 
quires a domain structure that is locally anisotropic and 
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FIG. 1. Off diagonal disorder and the anisotropy (a) A 

typical form of the employed off-diagonal disorder for L — 40. 
The off-diagonal component has a weakly disordered mod- 
ulation wavelength which is assumed to be close to the su- 
perconducting coherence length. In (b), the order parameter 
of the anisotropy as defined in the text, shows a strong a 
peak at 0.8t, and decays slowly at higher energies. The decay 
is associated with the fact that the domains become inter- 
nally homogeneous at very high energies, even though there 
is a remaining anisotropy. Given an experimental resolution, 
it is clear that decreasing the magnitude of A s makes the 
anisotropy observable at higher energies. 



highly similar to the observed STM patterns. The order 
parameter [l7| we use to identify the anisotropy, shows 
that by increasing the disorder amplitude, the anisotropy 
becomes noticeable at higher energies. Such anisotropy 
is not evident in a system with disorder that respects the 
d-wave symmetry. Finally, we show that disorder with ir- 
regular features, coming, for example, from an off-critical 
system, has no direct correspondence to the LDOS pat- 
terns, and at high energies there are no local anisotropic 
features, due to the complexity of the high energy spec- 
trum and induced bound-states [l8j]. 

To motivate the existence of modulated off-diagonal 
disorder, we consider a system which potentially can sup- 
port both s- and d- wave order parameters ^ s = |^ s |e 1( ^ s 
and = |^d|e 1( ^. The Ginsburg-Landau functional 
takes the form 

F[* d , *J = a d \* d \ 2 + ^|*,| 4 + K d \V* d \ 2 



+ a ad (9* 8 9 d + c.c.) + 6 fld |* fl | 2 |* d | 2 , (1) 

where a sd is the coupling between the two order param- 
eters and b sd > is a repulsive interaction that sup- 
presses all s-d mixing in the absence of the bilinear a sd 
term. Furthermore it is assumed that a d < a s < 0, 
b d , b s > 0, a sd . The energy of phase boundaries is 
F B = -J a |*i||*2|cos(^-^)-J ( ,|*J||*3|co S (^-^) ) 
where 1,2 label neighboring grains. In the limit where 
J s < Jd, b sd \ty d \ 2 > J s and a s <C a d , which is a nat- 
ural limit for the cuprates, minimizing the free energy 
favors generically the proliferation of phase domain walls 
of the s-order parameter where the b s d repulsion energy 
is suppressed [19], and the s-order parameter becomes 



zero on the boundary. Since the grain boundaries are de- 
fined within the coherence length scale, the off-diagonal 
order parameter A s should be modulated over distances 
larger than the coherence length. Additional motivations 
for the consideration of a modulated mixed SC order pa- 
rameter come from microscopically motivated mean-field 
states [l4j|, where the coexistence of modulated mixed 
order parameters is necessary. 

We consider a model of lattice itinerant fermions, in 
which the Hamiltonian takes the form, 

H F = -t ^(4r c j^ + h.c.) + ^(A ia c iT c i+c 4 + h.c.) 

(ij)cr i,a 

+ ^(A i (s) c iT c u + /i.c.) ) (2) 



where c\ G are fermionic operators, and Ai a =/?| A^e 2 ^ 1 are 
complex numbers for the d-wave SC order parameter de- 
fined now at the links (i,i+a)(a=unit vector along the x 
or y directions; f3=l (-1) for a along x (y)). The order 
parameter A is chosen to be modulated in space and 
take the values ±| A^ s ^ | in short - range domains [20 ]l(a)| 
For A( s ) — >• this model can be studied in perturba- 
tion theory. If we assume that A( s \x,y) = sin(</> g r), the 
spectral function at first-order takes the form, 

A(cj,q) = iGi T (k,k,-^) + iG u (k,k,cj) 



a j_ E l Sin Ok+q + </>k) 1 

^ sS * N 2^ 



,(3) 



which has poles at the same locations as the unperturbed 
system. This result shows that for weak disorder ampli- 
tudes, the low energy properties of the superconductor 
and the positions of the nodal points are unaffected, as 
it is a well known experimental fact for the cuprates [2l| . 

In the following, we go beyond perturbation theory 
and we study both low and high energy features of the 
model by solving it in a self-consistent mean-field man- 
ner [22[, using exact diagonalization on systems of size 
L x L. The relevant parameters to our calculation sat- 
isfy the hierarchy A^ < A^ < t, with t ~ 400raeV 
and A d ~ 60 — lOOmeF, which we believe to be consis- 
tent with the experiments on cuprates. Equation [2] can 
be diagonalized by using a Bogoliubov transformation. 
The corresponding Bogoliubov - de Gennes equations are 
solved iteratively, for fixed A^ s ^ fixed until a self consis- 
tent solution is found for the d-wave order parameter 




(4) 



The mean-field parameters are updated after each itera- 
tion using the equation 



A m = 9 ^2(u nji v^ i+a + Unj+sv*^) tanh ^ 
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FIG. 2. LDOS at Low energies. The typical form of the 
LDOS at a low energy uo = 0.3 is shown, for the pure 
5ao density wave model, and where a potential impurity is 
added, in order to compare with the results of Ref. |23j |. The 
amplitude of the Fourier transform of the patterns, shows 
a number of peaks that are not dispersing with the energy 
change, but their amplitude is weak and close to the experi- 
mentally resolved q2 wavevector [1]. The origin of the peaks 
is associated with the allowed scattering wavevectors for a d- 
wave superconductor in the presence of a density wave, similar 
with Ref. [23[ • In that case, the wavevector is (drzr/4, 0) and 
(0,±7r/4), here it is (±tt/5, ±3tt/5) and (±3tt/5, ±tt/5), very 
close to what has been labelled as qi. 



where the nearest neighbor interaction g < has to be 
attractive in order for the d-wave superconductor to be 
stable. The coupling g was chosen g = go = — 2M so 
that the self-consistent average d-wave gap A^ is 0.26t. 
The amplitude of the s-wave disorder is chosen as 0.05t, 
unless stated otherwise. 

The spectral function of the model (cf. Fig. |l(a)[ ) at 
low energies is similar to the unperturbed d-wave SC, 
where the low-energy quasiparticles extend along the un- 
perturbed fer mi- surface. In order to study the features 
of the LDOS at low energies in comparison to the well- 
studied unperturbed case 12|, 23[, we consider an impu- 
rity potential (a change in the chemical potential) at a 
single site with small strength U = O.Obt in order not 
to induce any Friedel oscillations. As shown in Fig. 
the LDOS contains similar peaks as the unperturbed 
system [l2[. In addition, there are more non- dispersing 
peaks at Q = (±7r/5, ±37r/5) and (±37r/5, ±7r/5), which 
vanish at larger energies, similar to the experimental ob- 
servations These peaks signify the emergence of the 
locally anisotropic density wave, which at low energies 
has very weak features. In the presence of spatial disor- 
der, like for example in the model of Fig. 1, these peaks 
are in close proximity to the denoted wavevectors Q. 
Since these peaks are static, they might be present in the 
region around the observed wave- vector qi [1]. The pres- 
ence of these disorder peaks could validate the present 
model. 

The weak disorder we described is adequate to induce 
local anisotropic features at high energies. As shown in 



Fig. 3(a), the LDOS shows strong anisotropic features 
and internal domain structure that resembles the ex- 
perimental observations. The large intensity near the 
domain walls on line-like structures and the similar- 
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FIG. 3. LDOS at high energies. The form of the LDOS at 
uj — 0.7t = 2.0 Ad is shown in (a). The pattern that emerges 
has strong anisotropic features, similar to the experimentally 
observed ones. Notice the bond-based structure and the local 
anisotropy that forms in this case. Such features are robust 
under weak disorder. In (b), the form of the LDOS when 
the disorder respects the d-wave symmetry, has no evident 
anisotropic features at uj = t. Only a natural modulation 
of the density of states is visible, imposed by the coupling 
modulation. 



ity to large parts of the observed STM maps, signify 
that our toy description captures important experimen- 
tal facts. Moreover, we define a local order parameter of 
the anisotropy [17|], applied on the LDOS F(r) 

Q(r) = {{dl - d 2 y )F{r)f + A{{d x d y )F{r)f (6) 

This order parameter, averaged spatially and over 
several disorder configurations, distinguishes between 
anisotropic fluctuations of the LDOS. As shown in 
the anisotropy increases monotonically until 



Fig. [1(b) , 

E t h — O.St and then decays slowly at higher energies. The 
energy threshold where the anisotropy is maximally vis- 
ible, is almost independent of the amplitude of the disor- 
der, but dependent on the d-wave gap scale (Eth ~ 2A^). 
Assuming that the experimental resolution allows for a 
low threshold on the identification of the order parameter 



(horizontal line in Fig. |l(b) ) indicates that the smaller 
the amplitude of the off-diagonal disorder, the higher 
the energy where the locally anisotropic features of the 
domains become visible. This observation signifies that 
the energy scale (~ 2A^) where the high-energy domains 
in STM experiments become visible can be much larger 
than the actual energy scale which leads to their forma- 
tion, which it is just 0.3A^. 

Another form of possible superconducting disorder is 
a pair-density wave (PDW) such as the one proposed in 
association to experiments in La2- a 3a :c Cu04 [24J and 
Bi2Sr2CaCu208+(5 [25[] . We investigated the behavior of 
LDOS in the presence of such a PDW, by setting A^ s ^ = 
and modulating the coupling g = go+5g, where Sg follows 
the pattern of Fig. |l(a)| an d has a maximum \Sg\ max = 
O.Obt. As seen in Fig. 3(b), there is no evidence of local 
anisotropy due to the presence of such disorder. 

Given that off-diagonal disorder generates high en- 
ergy bound states which distort the LDOS in complex 
ways [18], we studied the behavior of the LDOS in cases 
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FIG. 4. Stability of the anisotropy under strong spatial 
disorder. A typical off-critical (T = 1.4T C ) configuration of a 
conserved order parameter Ising model is shown in (a), which 
represents a typical form of strong spatial disorder. The form 
of the LDOS (only the upper left quarter of the configura- 
tion is shown) in (b) , at uj = t shows that the correspondence 
between the order parameter configuration and the form of 
the LDOS becomes very complicated. At very high energies 
(uj = t), there is no evident correspondence between the po- 
sitions of the domains and anisotropic features of the LDOS, 
due to the presence of a large collection of bound states gener- 
ated by the jerky features of the order parameter, as discussed 
in Ref. [11 



where the spatial form of the disorder has jerky features 
or where the amplitude has strong variations. When the 
defined domains have jerky features, as these would ap- 
pear, for example in a conserved-order parameter Ising 
model near its critical point, we find that the LDOS 
high energy features have complex characteristics with 
no local anisotropy at regions where large domains exist 
(cf. Fig. 4(a)|4(b)| ). However, amplitude (with no spa- 
tial jerkiness) variations of the off-diagonal disorder, as 
soon as they have zero global average, do not affect our 
qualitative conclusions. 

In conclusion we studied an explicit model of itin- 
erant fermions which provides a simple explanation of 
the experimental facts observed in recent STM experi- 
ments. This model contains off-diagonal disorder, which 
is known to be ubiquitous in systems with tendency to- 
wards orthorhombic distortions (like Y2_ x Ba cc Cu04and 
La2-xSr x Cu04), and in other less distorted systems, such 
as Bi2Sr2CaCu208+5and Ca2- x Na :c Cu02Cl2, it may oc- 
cur in the form of an inhomogeneous A s , modulated over 
distances of a few lattice sites. Off-diagonal disorder of 
this type does not distort the form of the LDOS at large 
distances and low energies, in contrast to potential im- 
purities [26[ . Furthermore its effect on the features of the 
LDOS, should not depend on the mechanism that drives 
them, which could be purely electronic as discussed in 
Ref. [3. 
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